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Tomograms are obtained as probability distributions and are used to reconstruct a quantum state from exper¬ 
imentally measured values. We study the evolution of tomograms for different quantum systems, both finite 
and infinite dimensional. In realistic experimental conditions, quantum states are exposed to the ambient envi¬ 
ronment and hence subject to effects like decoherence and dissipation, which are dealt with here, consistently, 
using the formalism of open quantum systems. This is extremely relevant from the perspective of experimental 
implementation and issues related to state reconstruction in quantum computation and communication. These 
considerations are also expected to affect the quasiprobability distribution obtained from experimentally gener¬ 
ated tomograms and nonclassicality observed from them. 
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I. INTRODUCTION 


A quantum state can be characterized by a number of probability and quasiprobability distribution functions (1. The 
quasiprobability distributions are not true probability distributions as most of them can have nonpositive values. Interestingly, 
this nonpositivity can be viewed as a signature of nonclassicality or quantumness. Specifically the negative values of the Wigner 
H and P i,il functions serve as witnesses of nonclassicality. Further, zeros of Q function @] also serve as witness of nonclas¬ 
sicality. As there does not exist any straight forward prescription for direct measurement of these quasiprobability distributions, 
several efforts have been made to construct measurable probability distributions that can be used to uniquely construct either all 
or some of these quasiprobability distributions. Such measurable probability distributions are referred to as tomograms 
In other words, the tomogram is a scheme for measuring a quantum state by using a representation in one to one correspon¬ 
dence with the true probability distribution rather than with a quasidistribution EE A relationship between a tomogram and a 
quasidistribution function, such as the Wigner function, can be established for both continuous and discrete systems El El. 
Specifically, in Ref. iflUl it was shown that quasiprobability distributions (P, Q, and Wigner functions) can be uniquely deter¬ 
mined in terms of probability distributions for the rotated quadrature phase which can be viewed as an optical tomogram of the 
state. Similarly, in Ref. 0 it was shown that for finite dimensional phase states, discrete Wigner functions and tomograms are 
connected by a discretization of the continuous variable Radon transformation and was referred to as the Plato transformation. 

In the recent past, a few successful attempts have been made to measure Wigner function directly in experiments lfl3lll4l] . but 
the methods adopted are state specific. The same limitation is also valid for the theoretical proposals lfl5ll for the measurement 
of Wigner function. Further, optical homodyne tomography has been employed for the experimental measurement of the Wigner 
functions of vacuum and squeezed states in lEbiEl . while distributions corresponding to Pegg-Barnett and Susskind-Glogower 
phase operators were also obtained in 0 . An experimental measurement of the P, Q and Wigner quantum phase distributions 
for the squeezed vacuum state has been reported in flUl . Precision of homodyne tomography technique was compared with 
conventional detection techniques in 0 . A number of alternative methods of tomography have also been proposed (2Cj-[22j] , and 
exploited to obtain phase distributions like Wigner and Q functions [23]. In 124 1 continuous variable quantum state tomography 
was reviewed from the perspective of quantum information. In brief, there does not exist any general prescription for direct 
experimental measurement of the Wigner function and other quasidistribution functions. In practice, to detect the nonclassicality 
in a system the Wigner function is obtained either by photon counting or from experimentally measured tomograms f IS]. Thus, 
tomograms are very important for the identification of nonclassical character of a physical system. In another line of studies, 
simulation of quantum systems have been performed using tomography. For example, tomograms were used for simulation of 
tunneling |25-127|] and multimode quantum states lf28ll . Attempts have also been made to understand the tomogram via path 
integrals 1291 1301. 

Furthermore, how to reconstruct a quantum state from experimentally measured values is of prime interest for both quantum 
computation [130 and communication Blil. Specifically, in Ref. 0 it is strongly established that tomography and spectroscopy 
can be interpreted as dual forms of quantum computation, and in Ref. 0, quantum teleportation was experimentally performed 
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over a distance of 143 km and the quality of teleportation was verified with the help of quantum process tomography (QPT) of 
quantum teleportation without feed-forward. Here it would be apt to note that OPT is an aspect of quantum state tomography in 
which a quantum process is obtained as a trace preserving positive linear map [ 32 L In the recent past, quantum process tomogra¬ 
phy has been discussed from the perspective of open quantum system effects 133143511 . A novel method of complete experimental 
characterization of quantum optical processes was introduced in lf36tl . It was further developed in (37, 38] and extended to 
characterization of iV-modes in S. In El, QPT was applied to the characterization of optical memory based on electro- 
magnetically induced transparency while 14111 and |42] were devoted to QPT of the electromagnetic field and conditional state 
engineering, respectively. Quantum state tomography has its applications in quantum cryptography as well El - Specifically, in 
Ref. El an interesting protocol of quantum cryptography was proposed in which eavesdropping in the quantum channel was 
checked by requiring consistency of outcome of the tomography with the unbiased noise situation. Keeping these facts in mind, 
we aim to construct tomograms for a number of physical systems of practical relevance (mostly having applications in quantum 
computation and communication) and investigate the effects of various types of noise on them. 

From the experimental perspective, a quantum state always interacts with its surroundings. Hence, the evolution of the cor¬ 
responding tomogram after taking into account the interaction of the quantum state with its environment should be considered. 
This can be achieved with the open quantum system formalism EH-El- Specifically, both purely dephasing (QND) Ei and 
dissipative El open quantum system effects have been studied here. Interestingly, both these effects have also been experimen¬ 
tally realized in the recent past |[50i . l5ltl . In Ref. ['52}], a systematic study of quasidistribution functions was made for a host of 
interesting states under general open system evolutions. 

Here, we set ourselves the task of obtaining the tomograms for various finite and infinite dimensional quantum systems in 
different open quantum system scenarios. For finite dimensional spin states, the tomogram is the distribution function of the 
projections of the spin on an arbitrary axis, characterized by Euler angles, and can be obtained from the diagonal elements of 
the rotated density matrix, while for continuous variable systems, such as the radiation field, the analog would be the homodyne 
probability. It follows from general group theoretical arguments that, making use of unitary irreducible square integrable rep¬ 
resentation of the tomographic group under consideration, a unified tomographic prescription can be developed for both finite 
dimensional and continuous variable systems |53]. Tomograms for spin states have been developed both as projections on an 
arbitrary axis [54] as well as by using_a discrete variable analog of symplectic tomography lEsi. Tomograms of optical systems 
have been well studied in the past lla. IllL[24.I56L 157 1. In Ref. [58] quantum state tomography was used to determine the degree 
of non-Markovianity in an open system. Further, thermal noise is used in tomography (for reconstruction of photon number 
distributions) as a probe |59j|]. 

The paper is organized as follows. In Section [TI] tomograms of single spin--} (qubit) atomic coherent state under purely 
dephasing (QND) and dissipative evolution are obtained. Further, the tomogram of two spin-} (qubit) quantum state is studied 
in Section Hill under the influence of a vacuum bath. This is followed by a tomogram for a general spin-1 pure state in Section 
ITVl The tomograms of finite dimensional number-phase states under open quantum system evolution are discussed in SectionlVl 
This is illustrated by a specific example of a three-level quantum (qutrit) system evolving under a spontaneous emission channel. 
In Section IVH we discuss the tomogram of an infinite dimensional system, the ubiquitous dissipative harmonic oscillator. We 
conclude in Section lVTII 


II. TOMOGRAMS OF SINGLE SPIN-} STATES 

In this section, we study the tomograms for single spin-} (qubit) atomic coherent state evolving under two general noise 
models, i.e., pure dephasing (QND) and dissipative squeezed generalized amplitude damping (SGAD) evolution, incorporating 
the effects of dissipation, decoherence and bath squeezing. 


A. QND Evolution 

The master equation of a quantum state under QND evolution EH] is 

Pnm (t) = -I ( E n ~ E m ) + if] (t ) (E% - E^) - ( E n - E m f 7 (t) p s nm {t) , 
where E n ’s are the eigenvalues of the system Hamiltonian in the system eigenbasis | n), 

V(t) = sin ( Wfct )’ 


( 1 ) 


and 


7 (t) = lY.W-z coth (| {e iuikt - 1) cosh (r k ) + ( e~ iUkt - l) sinh (r k ) e 2i$fc 

k k ' ' 
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with f3 = -fT^r- Here, kn is the Boltzmann constant, 77 and <1>are the bath squeezing parameters and is the system-bath 
coupling coefficient. The initial density matrix for the atomic coherent state is given by 


p s (0) = \a, /3|, 


( 2 ) 


where the atomic coherent state is given by 


K/ 3 ) = E 

m——j 


2 j 

j + TO 


1/2 


sin (f y +m cos (f y m I j,m)e~ i ^ +m ^. 


The different elements of the density matrix in Eq. (0 at time t under QND evolution becomes 

s (j\ _ -iu](m-n)t i(hui ) 2 (m 2 -n 2 }r)(t) -(hw ) 2 (m—n) 2 ^(t) n s /o') 

rim,in \ L J c c c rim,in \ u / i 


(3) 


(4) 


with Pj m j n (0) = (j, m\p s (0) \j, n). Considering the initial state of the system as atomic coherent state, i.e., using Eq. (0. 
different elements of the density matrix in Eq. 0 at time t = 0 are 


Pjm,jn (0) 


2 j 

j + to 


1/2 


2 j 

j + n 


1/2 


,i(n—m)l3 l a\ 2 l+ m + n 


Sill 


(tr 


cos 


(f) 


‘Ij—m—n 


(5) 


Using Eq. 0 as the initial density matrix elements in Eq. 0 , we can write all the elements of the density matrix at time t as 

I/ 2 / r, ■ \ t/ 2 


Pjm,jn (^) ^ 


2 j 

j + m 


2 j 

j + n 


0 —iui(m—n)t„i(hui) 2 (m 2 —n 2 ^)r)(t) 


( 6 ) 


x e -(M a (m-n}» 7 (t) sin (a)«+"*+" CQS (| fj-m-n 


To obtain a tomogram of a spin-2 atomic coherent state under QND evolution, we can express the density matrix in terms of 
Wigner-Dicke states as 


P U) = P U) (t) = Pm,m’\j> m ){j> m '\- ( 7 ) 

m,m'=—j 

The different elements of this density matrix p^ m , = (m|p( J ) | TO ') ca n be obtained using Eq. 0 , with 77 i, n —> to, m! , for 
j = 2, to, m! = ±2. Subsequently, the density matrix is obtained as 


P 


( 1 / 2 ) _ 


Sln2 (f) 2 

\ e iujt e -(hu,) 7(t) s i llae */3 


l e -iu,t e -(hu 1 ) 2 1 (t) sinae -i/3 
2 


COS- (f) 
1/2 


( 8 ) 


We can easily check that the trace of the density matrix (p*- 1 / 2 -*) is one, i.e., PnJ.m = 1. Further, the tomogram of 

m=-1/2 

this state can be expressed as 


where D 


lo (toi, a, f 3 , 'y) = £ E (“-M) » 

' / m=—jm'=—j v ' 

{j \r n , (a, is the Wigner //-function 


(9) 




,—im 7 


( 10 ) 


and the notation used here is consistent with that in Ref. a. Here, a, /?, and 7 are Euler angles = <p , 9, and ip, with 
<p,ip € [0, 27t], and 0 £ [0,7r], and 


d^m' (?) 


(j+m)\(j-m)\ 


1/2 / m+m' 

f ■ b\ 

(cos f J 

(sm f J 


( 11 ) 
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Figure 1: (Color online) The variation of the tomogram with time (t ) for single spin-1 atomic coherent state in the presence of QND noise with 
bath parameters 70 = 0 . 1 , co c = 100, squeezing parameters r = 0, a = 0, and u> = 1.0 and a = ft = jt, /3 = 7 = tt, in the units 

of h = kB = 1. The smooth (blue) line, dashed (red) line and dot-dashed (magenta) line correspond to the tomogram with time for different 
temperatures T = 0, 1 and 2, respectively. 


where Pi a,b ' 1 ( x ) are Jacobi polynomials. A tomogram is the spin projection onto an arbitrary, rotated, axis. The physical 
significance of the ZZ-function is its connection to the process of rotation and can be illustrated by 


|j,mi) = D { ^ iin ,(a,p, j) , 

{j,m' 2 \Rl (a,£,7) (a, / 3 ,7J . 

Here, R (a, /3, 7 ^ stands for the operation of rotation about an axis whose orientation is specified by 5, 3- and 7 . Using the 
different values of m and m’, we can obtain various Wigner ZZ-functions as 

fl | 1/2 ! (“.ftl) = - sin Tf) e “ 5 ( s -t), 

= cos (“+=(), 

(a,/3,7^ = cos(|je5(“+7), 

ZZifi ( 5 ,^, 7 ) = sin e5( Q “7). 

Using the hrst two relations of Eq. ( IT3l > and Eq. <[8]», the first component of the tomogram can be obtained from Eq. ([9]) as 


w a, /3, 7 ^ = wi = cos 2 — cos/3cos 2 (^) — | sin/3sinacos (wf + /3 + 7 ) e ( K “) 2 7 ( 4 ). ( 14 ) 

From Eq. (fTTb . it can be inferred that the tomogram is free from Euler angle a, and consequently is a function of 3 and 7 only, 
or /(/3, 7 ). It is worth mentioning here that 7 and 7 (t) are two different parameters, the former being an Euler angle while the 
latter is responsible for decoherence. The variation of the tomogram is given in Fig. [T]with time, for the different temperatures. 
For the second component of the tomogram with m\ = —7, using last two relations of Eq. (fT3l> and substituting Eq. ® in Eq. 
©. we obtain 

w (—^, 5 ,/Sjt) = W 2 = cos 2 (|) -cos/3cos 2 (f) +isin^sinacos(u;Z + ^ + 7 )e-^) 2 ^). (15) 


We can check the validity of the tomogram obtained by verifying that = 07 + UJ 2 = 1- Interestingly, we can see that 
the knowledge of one of the components of the tomogram is enough to reconstruct the whole state. Keeping this in mind, we 
have only shown the variation of 07 in Fig. [T|as 0 J 2 = 1 — oji- 

In Fig. [I] we can easily see the expected behavior of tomogram with increase in temperature for zero bath squeezing. Specif¬ 
ically, with increase in temperature, the tomogram tends to randomize more quickly towards probability 1/2. Fig. [2] further 
establishes the effect of the environment on the tomogram. Particularly, Fig. |2]b brings out the oscillatory nature of tomogram 
with time while temperature tends to randomize it. Similarly, Figs. [2] a and c show the dependence of the tomogram on Euler 
angles and the atomic coherent state parameters, respectively. 
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Figure 2: (Color online) The dependence of the tomogram on various parameters is depicted for single spin-| atomic coherent state in the 
presence of QND noise with bath parameters 70 = 0.1, u c = 100, squeezing parameter a = 0, and w = 1.0 in the units of h = ks = 1. In 
(a) the tomogram is shown as a function of f3 and 7 with a = ^, /3 = ^, and r = t = T =1; while (b) exhibits the variation of the tomogram 
with time and temperature for r = 0 and a= : |,/3= : |,/3=^,7=j. Finally, the last plot (c) shows the dependence of the tomogram on 
the atomic coherent state parameters a and /? with /3 = £, 7 = f at time t = 1 for bath squeezing parameter r = 1 at T = 1. 


B. Dissipative SGAD channel 


Master equation for the dissipative evolution of a given state in the squeezed generalized amplitude damping (SGAD) channel 
is given by [49fl 

j- t p s {t) = -if[a z ,p s (t)\+^(N+l){a-p s (t)cj + -\a + a-p s (t)~\p s {t)a + a-} 

+ {(J+p 8 (t) <7— - \<7-<j + p 3 (t) - \p s (t) cr_CT + } 70 N - 70 Ma + p s ( t ) a+ - 7 0 M*<j-p s (0) ct_. 

The density matrix for a quantum state under a dissipative SGAD channel at time t, can be obtained, from the above equation, as 


p s (*) = \p s (0) /+ + \<?zP s (0) <J z f— - JP S (0) cr z g_ - \a z p s (0) g+ - 7o^^e V 

x {Mcr + p s (0)a+ + M*a-p s (0)a-} + {^v<?-P s (0)cr+ + ^a+p s (0) cr_| , 


(17) 


where 


/ ± = { 


= U + ‘ ± 2 


cosh (a't) e V } , g± = (l - e ^ ) ± ^ sinh [a't) e Vj , 7+ = 7o (AT + 1), 7 _ = 


7 o 7 V, 7 ^ = 7 + + 7 _, 7 = 7 + - 7 - = 70 , a' = |M | 2 - w 2 ; 


and 


a+ = |1)(0|, a- = 10)(11, 

(7 z = 

= | 1 )( 1 | - | 0 )< 0 | 

= Ie)(e|-|p)( 5 |. 


Also, 


<?z\g) = -I g), Oz\e) = |e); 

v+\g) = |e), er+|e) = 0; 

cr-jp) = 0, cr_|e) = |p). 


Here, 70 is the spontaneous emission rate, M = —1 [2N th + 1} exp (i(f>) sinh (2 r), and N = N th {cosh 2 (r) + sinh 2 (r)} + 
sinh 2 (r) , where Nth = 1 / {exp (hui/ksT) — 1} being the Planck distribution, and r and the bath squeezing angle ((/>) are 
the bath squeezing parameters. The initial state, as for the tomogram of a quantum state under QND evolution, is the atomic 
coherent state given in Eq. (J2]». Using Eq. ( fTTI) . the density matrix can be written as 


P s C t ) 


aip s (*m> <iip s w 1 -1)' 

HIp s WI5> 


(18) 
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where the various terms are 

(W toll) 

(■ \\p s (f)l-!> 
{-\\P S Wl3> 

(~\\p s (t)l-l) 


sm 




{cosh (a't) - £ sinh (a't)} e~V - ^ sinh (a't) e# 

: sm a {cosh (a't) + ^ sinh (a't )} e I/3 — 7 ° < ^f sinh (a't) e~ 1 ^ 
cos 2 (f) e-^ + ^ (l - e" 7 '*) , 


1/2 

and the density matrix can be seen to be normalized as Y (m\p s (t) \m) = 1. 

m—- 1/2 

The tomogram of a state evolving in a dissipative SGAD channel, in analogy to the QND case, can be obtained by expanding 
the density matrix in the basis of the Wigner-Dicke states, as in Eq. CD- Using Eq. Q- the first two relations of Eq. (fl3l > and Eq. 
the first component of the tomogram is 



sin 2 | cos 

±sin/ 3 {e" i7 


2 (f) 

\ sinae _l/3 


c.c.}. 


+ 

e~ 


7+ 


0 

{cosh (a't) — 


+ cos2 (2) { sin2 (f) e (i - e 7/3 ‘) } 

^sinh(a'f)} - sin a sinh (a't) 


(19) 

Again, we can check the validity of the analytic expression of the tomogram in the absence of open system effects, i.e., by 
considering 70 = 7 = 0, 7+ = 7_ = 0 = 7^, which leads to a' = ico, we have 

= cos 2 — cos^cos 2 (^) — \ sin/ 3 sinacos (ojt + /3 + 7), ( 20 ) 

which is identical to the QND case, i.e., Eq. (f]~4l >. with 7 (t) = 0. Similarly, using Eq. ©. the last two relations of Eq. (IT3l) . and 
Eq. m, we obtain the second component as 


<*;(-§, a, £,7)=w 2 = cos 2 (cos 2 (f) e + 2+ (l - e 7 ' 3 ^ | + sin 2 (§) (sin 2 (f) e ^ (l - e 7/3 *( j 

+ }sin/ 3 (e - * 7 } sinae - * /3 e _2 2- {cosh (a't) — ^ sinh (a't)} — 7^- since sinh (a't) 

+ c.c.}. 

( 21 ) 

Similar to the first tomogram of the dissipative SGAD channel, we can check the solution in the absence of the open system 
effects which leads to a' = ioj. This can be seen to be the same as the corresponding QND case, i.e., Eq. 03, with 7 (f) = 0 


(— i,a,/3,7^ = cos 2 — cos/ 3 sin 2 (^) + i sin/ 3 sinacos (wf +/3 + 7). 


( 22 ) 


We can also check the validity of the tomogram as in the QND case by Y = w i + w 2 = 1. Hence, as before, one component 
of the tomogram would be enough to recover all the information. This is why, in the plots, we only show the first component of 
the tomogram. The other component can be easily obtained from it. 

The variation of tomogram with different parameters is shown in Figs. [3] and [4] Fig. [3] exhibits the randomization of the 
tomogram with increase in temperature. This fact can be observed in the smooth (blue) and dashed (red) lines. However, an 
interesting behavior is observed here with respect to bath squeezing. It can be seen that it takes relatively longer to randomize 
the tomogram in presence of squeezing than in its absence, temperature remaining same, as illustrated by a comparison between 
the dot-dashed (magenta) and dashed (red) line. This fact, in turn, establishes that squeezing is a useful quantum resource. This 
behavior is further elaborated in Fig. [4] where the effect of bath squeezing can be observed and is consistent with the quadrature 
behavior of squeezing. This beneficial effect of squeezing is not observed for evolution under QND channel. From the present 
analysis, it could be envisaged that a tomographic connection could be established between the state, under consideration and 
the generic open system interaction evolving it. 


III. TOMOGRAM OF TWO SPIN-± (QUBIT) STATES 


Various two-qubit tomography schemes have been proposed in the recentpast 061 1465 0. Specifically, the tomogram for two 
spin-} (qubit) states can be obtained using the star product scheme 1161 L 1621. In (640, two-qubit states were analyzed from 
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Figure 3: (Color online) The tomogram varying with time ( t ) is shown for a single spin-| atomic coherent state in the presence of the SGAD 
noise for bath squeezing angle <j> = 7 r in the units of ft = fcs = 1, with uj = 1.0, 70 = 0.25, and a = j, 0 = 0 = 7 = j. 

The smooth (blue) line, dashed (red) line and dot-dashed (magenta) line correspond to the tomogram with time for different temperatures and 
squeezing parameters T = 1, 10 and 10, and r = 0, 0 and 1, respectively. 



60 -1 


Figure 4: (Color online) The tomogram for single spin -7 atomic coherent state in the presence of SGAD noise is shown as a function of 
squeezing parameters r and <f> with a = f3 = ^,/3 = f ,7 = 7 , and uj = 1.0, 70 = 0.25, in the units of ft = kn = 1 for T = 1 at time 
t = 1 . 


the perspective of tomographic causal analysis, while in |65}l, an interesting connection between tomographic construction of 
two-qubit states to aspects of quantum correlations such as discord and measurement induced disturbance was developed. 

For a two qubit state p one can obtain the tomogram as 


uj (mi, m2) = Tr \p{Qi (mi) <g> Q 2 (m 2 )}\ , (23) 

where Qi (mi) = Uj \ m;) (mi\ Ui, and m; = ±|, while the unitary matrices Ui are 


cos Tr exp ■ 

i(a«+7«) 1 

sin 4r exp 

f i(ai- 74 ) ) 

2 1 

L 2 i 

• Bo 

i(«i— fi) 

2 

1 fa / 


— sm ^ exp 

> cos exp < 

2 J J 


for i £ {1,2}. Hence, the tomogram of the two qubit state can be written as the diagonal elements of p, where p = 

(Ui 0 U 2 ) p (Ui g> U 2 )^. 
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A. Tomogram of two qubits under dissipative evolution in a vacuum bath 


Now, we construct the tomogram of a two qubit state in a vacuum bath under dissipative evolution, as discussed in Ref. 

The initial state of the system is considered with one qubit in the excited state \e\) and the other in the ground state \g 2 ), i.e., 
I e i) |ff 2 )- The reduced density matrix of the system of interest, here the two qubits, is 


Pee (t) Pes ij') Pea (t) Peg (t) 
PesW Pss(t) Psa(t) Psg(t) 
Pea W Psa W Paa (t) p ag ( t ) 
Plg{t) P*sg(t) P*ag(t) Pgg (t) 


(24) 


P(t ) = 

where the analytic form of all the elements of the density matrix is given in Appendix 1. 

The tomogram can be thought of as a tomographic-probability vector ui = [tui, u> 2 , W 3 , W 4 ] 7 " (here T corresponds to transpose 
of the vector), where each component can be expressed analytically as 

UJ 1 (t) = \ 4p ee cos 2 If cos 2 If + 4 p gg sin 2 \ sin 2 If + ( p aa + p ss ) (l - cos /3i cos /? 2 ) 


(Paa ~ Pss) sin/3i sin/3 2 cos (71 - 72 ) + jp SQ (cos/3i - cos (3 2 - i sin^i sin/3 2 sin (71 - 7 2 )j 


+ y/2 


-Pea + Pes) cos 2 ^ + ( p ag + p sg ) sin 2 ^ j sin/3i exp ( 271 ) + sin/3 2 exp ( 772 ) 


(25) 


X ((Pea + Pes) cos 2 ^ - ( p ag - p S g) sin 2 ) + exp (771 + *72) p eg sin / 3 i sin / 3 2 + c.c.j 

W2 (f) = \ 4 p ee cos 2 \ sin 2 If + 4 p gg sin 2 If cos 2 If + ( p aa + p ss ) (l + cos ^1 cos / 3 2 ) 

+ (Paa - Pss) sin / 3 i sin / 3 2 cos (71 -72 ) + jp sa (cos ft + cos ft + i sin ft sin ft sin (7! - 72 

+ ((p ag + p sg )cos 2 If - (p eQ - p es )sin 2 |f j sinft exp(i7i) 

+ sinft exp (772) (- (p ea + Pes) cos 2 If + ( p ag - p sg ) sin 2 If j - exp (471 + *72) p eg sinft sinft + c.c.j 

w 3 (f) = 7 4 p ee sin 2 I s - cos 2 If + 4 p gg cos 2 If sin 2 If + (p aa + p ss ) (l + cos ft cos ft) 

+ (Paa - Pss) sin ft sin ft cos (71 - 72) + j-p sa (cos ft + cos ft - i sin ft sin ft sin (77 - 72) ) 

+ a /2 (- (Pas + Psg) Sin 2 If + (p ea - p es ) cos 2 If) sin ft exp (171) 

+ sin ft exp (172) ((p eo + p es ) sin 2 |f - (p ag - p sg ) cos 2 If ) - exp (*7! + 172) p eg sin ft sin ft 


(26) 


c.c.j 


(27) 


and 

w4 (i) 


4p ee sin 2 |f sin 2 |f 


+ 4p gg COS 2 \ COS 2 If + (p aa + p ss ) (l - COS ft COS ft ) 

(Paa - Pss) sinft sinft COS (71 - 72) + j psa (cos ft - cos ft + i sinft sin ft sin (71 - 72) ) 

+ v/2 (- (p ag + p sg ) cos 2 If + (p ea - p es ) sin 2 |f) sinft exp (ryi) 

+ sinft exp (ry 2 ) (- (p eQ + p es ) sin 2 \ + (p ag - p sg ) cos 2 If) - exp (*77 + *72) p eg sinft sinft 


c.c.j 


(28) 

Here, pi 3 are the elements of the matrix in Eq. (l24l > and are given in Appendix 1. For simplicity of notations the time dependence 
in the arguments of matrix elements is omitted. Similar to the tomograms for single spin-4 states the tomogram obtained here is 
also free from a. 

As in the cases of single qubit tomograms, we can again verify that the tomogram obtained here satisfies the condition 
Y, 0 Ji = p ee + p gg + paa + Pss, which is the trace of the density matrix given in Eq. (l24l) . and hence equal to one. 

For the case of identical qubits considered here, we take the wave-vector and mean frequency to be ft = ujo = 1, the 
spontaneous emission rate ft = 0.05 and p ■ fij = 0. Here, ji is the unit vector along the atomic transition dipole moment and 
fa is the inter-atomic distance. Further, the initial state of the system is taken to be p ee (0) = p oa (0) = p es (0) = p ea (0) = 
(0) - (0) = (0) - 0 . and ( 0 ) = „„„ ( 0 ) = (0) = 0.5. 

The variation of all four components of the tomogram is shown with different parameters in Figs. 0 and [ 6 ] In Fig. 0 
large oscillations can be observed for small interqubit spacing, which is consistent with the earlier observations in a plethora of 
scenario (5jE[66l[67[|. Fig. 0further demonstrates similar behavior for small interqubit spacing. For small interqubit spacing the 
ambient environment opens up a channel between the qubits resulting in enhancement of oscillations. 
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Figure 5: (Color online) Various components of the tomogram changing with time are shown in (a)-(d) for the two-qubit state, in the presence 
of vacuum bath, with /3i = 71 = ^, 72 = j and the inter-qubit spacing m = 0.05 (2.0) corresponding to smooth blue (red 

dashed) line. 


IV. TOMOGRAM OF SINGLE SPIN-1 STATE 


The tomograms for finite spin states have been considered, among others, in Refs, 
theme of this work, we take up an arbitrary spin-1 state 


where N = 


\/|a | 2 + |&| 2 -t-|c|' 


= N 

is the normalization factor. The corresponding density matrix is 


pW = |/V | 5 


|a| ab* ac* 
a*b \bf be* 
a*c b*c |c| 2 


. In continuation with the 


(29) 


(30) 


Here, we restrict ourselves to obtaining the tomogram for the state ( l30l ). without considering open system effects. Using Eqs. 
dOt- llTll . all the Wigner /9-functions for the tomogram can be calculated as before. Using Eqs. © and ( |30| ). ui ^ 1,5, /3. can 
be written as 




= wt = | Nf 


{4- (l +cos/3) 


+ 4 sin 2 /3 - 


kl ; 


{[-^)smp{l + cosp)+^f 1 sm 2 p 


(l-cos/3) | 


bc*e'2. 

2^2 


Sill 


/3 (l — cos /3) + c.c. |J 


(31) 


Similarly, w ( 0 , a , /3, 7 ) 


can be obtained as 
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4 


r 12 


Figure 6 : (Color online) (a)-(d) depict the tomogram for the two qubit state, interacting with a vacuum bath, as a function of the inter-qubit 
spacing at t = 1 (smooth blue line), and t = 5 (red dashed line). For all the plots /3i = /?2 = 71 = 72 = j- 


w(o,a, /3,7^j = w 0 = |1V| 2 sin 2 (3 + |6| 2 cos 2 /3 + sin 2 /?) 


2 

+ {s|^Zsin2^-s^f^sin 2 ^-^^sin2^+c.c.} 


(32) 


Also, w (—1; a, /3, 7^ is 


(-l,a,/3,7) 


= w-t = |(V| 

_i_ / _ 

+ I 2V2 


^j- (l - cos/3) + sin 2 /3 + (l + cos/3^ 


sin/3 (l — cos (3^ + ac sin 2 f3 + bc *^ sin/3 (l + cos/3^ + c.c.j 


(33) 


Interestingly, the tomogram obtained for a general spin-1 quantum state is also free from a as for spin-4 cases discussed 
above. Further, it can be checked here that the tomogram satisfies the condition u>i + ojm + w_ 1 = 1. For a = l. b = 0 = c, 
the tomogram, obtained here, is seen to be consistent with the results reported earlier | 683 > 


w (l,5, /3 ,t) 
w (o,5, /3, 7 ) 
w (-1,5, 


(l+COS/3) 2 
4 5 

(l—cos 2 /3) 

2 ’ 

(l — COS y0) 2 

4 • 


(34) 


The variation of all three components of the tomogram with Euler angles is given in Fig. [7] The peaks in one component have 
corresponding valleys in other components of the tomogram, which are manifestations of normalization of the tomogram to one. 
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Figure 7: (Color online) The dependence of all three components of tomogram for single spin-1 state with Euler angles for a = b = c = y=. 


V. TOMOGRAM OF A FINITE DIMENSIONAL STATE 

In this section, we discuss tomography of finite dimensional states. The tomogram of a finite dimensional state can be defined 
as m 


d -1 

UJ (m, t, q) =y] W ( tm - qx, qm + 1\) , (35) 

x=o 

where W (x, m) is the discrete Wigner function with integers t and q, depicting complementarity between number m (denoting 
angular momentum) and phase x °f finite dimensional states. Here, the discrete Wigner function is expressed as 


d -1 


W (x, m) = i ^2 exp ( ym0 ) (x ~ ©I P \x + ©) • 


e=o 


(36) 


It would be apt here to mention that the phase states |x) are periodic, such that |x + d) = |x). A general d dimensional density 
matrix can be written in Weyl operator basis as IH 

^ d—1 

P = ^ 1 + ^ bnmUnm ( 37 ) 

n,m—0 

d—1 

with boo = 0 and U nm = ex P T an ) l a ) ( a + m l- m ake use of the periodicity of phase states in computations 

a—0 

involving U nm . 

Up to now the treatment is applicable to any generic finite dimensional system. Here, for concreteness, we concentrate on an 
important finite dimensional system, viz. a qutrit with d = 3. We study the effect of spontaneous emission (SE) channel 0] on 
the qutrit. SE is a dissipative process which can be modeled by the following Kraus operators 


K 0 


'10 O' 

771 t 

0 e~+ 0 

770 1 

0 0 e~+ 


K x 


0 y/1 - e- 7 ? 1 * 0 
0 0 0 

0 0 0 


K 2 


0 0 v/1 - e 

0 0 0 

0 0 0 


(38) 
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Figure 8: (Color online) The variation of different components of tomogram for a qutrit state with time in the presence of a spontaneous 
emission channel with Einstein coefficients 771 = 2 and 772 = 4 . The smooth (blue), dashed (red) and dot-dashed (magenta) lines correspond 
to u>m with in = 0, 1, and 2, respectively. 


where 771 and 772 are two Einstein coefficients which control the population of the excited states. Thus, we can write the density 
matrix of an arbitrary three dimensional state at time t evolving under the spontaneous emission channel as 


Using this we can obtain 


p M = E (*) Ij> < fc i =E K *p (°) K l 


j,k =0 


i=0 


W (x, m, t ) 


1 

3 


2 

E A x-e, x +e (*) exp 

0=0 



( 39 ) 


( 40 ) 


As mentioned above x — © and x + 0 are mod d operations. Here, we have considered an initial density matrix given by Eq. 
( | 37 | > with 601 = bio = \ and 612 = 621 = while the remaining coefficients can be obtained from these values. Hence, the 
obtained tomogram with t = 0 and q = 1 has three components as 


uj ( 0 , 0 , 1 ) = w 0 
w (1,0,1) = u 1 
w (2,0,1) = u> 2 


30 

J_ 

60 

_1_ 

60 


10 + 7 i + e 2 ^ + e 5 ; 
20 - (e-^ - 13e~i^ 1+r ' 2 '> t , 

20 - 13 (e-^ + e-^ + lie ~h(m+V 2 )t 


( 41 ) 


2 

It can be easily seen here that the tomogram obtained is normalized as ^ co m = 1 in Eq. fiTl) . 

m —0 

For specific values of parameters of the spontaneous emission channel, i.e., Einstein coefficients, the evolution of tomogram 
for the qutrit state is shown in Fig. [8] The tomogram shows that the noisy channel tends to randomize all the components of the 
tomogram to one-third. We have already noted that once a tomogram is obtained for a finite dimensional system, it is possible 
to transform it to obtain the Wigner function for the system, and vice verse, but in an experiment we obtain tomograms. A lot 
of work has been devoted to the study of Wigner functions for finite dimensional systems BIziQ I- Some efforts have also 
been made to study tomograms for finite dimensional coherent states lF 7 l IbL p 74 h . However, to the best of our knowledge, no such 
efforts had yet been made to study evolution of tomograms in noisy environment. 


VI. OPTICAL TOMOGRAM FOR A DISSIPATIVE HARMONIC OSCILLATOR 

In the end, we come to the to mog ram of an infinite dimensional system, the harmonic oscillator. This is typical of a plethora of 
oscillatory and optical systems [ 44 , 75 ]. In Ref. E 3 ], the quantum mechanics of the damped harmonic oscillator was examined, 
from the perspective of a classical description of quantum mechanics & Use was made of the generating function method, 
resulting in the avoidance of the need to evaluate the Wigner function as an intermediary step for obtaining the tomogram. 
Further, in fzi the density matrix, state tomogram and Wigner function of a parametric oscillator were studied. 
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(a) (b) (c) 


Figure 9: (Color online) The effect of interaction of the optical tomogram with its environment is shown as a function of X and 9 for 
N = 5, r = 1 for the initial coherent state parameter /3 = 2 at time t = 0, 1, 10, respectively. 


Here, we construct the tomogram of the dissipative harmonic oscillator evolving under a Lindbladian evolution, in a phase 
sensitive reservoir 0. It would be pertinent to mention that tomographic reconstruction of Gaussian states evolving under a 
Markovian evolution has also been considered in 1341. The dissipative harmonic oscillator can be described by the Hamiltonian 

H = H s + H r + H S r , (42) 


where the system Hamiltonian H s of a harmonic oscillator is described as 

TT P 2 , 1 2 2 

Hs= 2^ + 2 mWX ' 

while the reservoir Hamiltonian Hr is given by 


h r = J2 


2 rrij 




with the system-reservoir interaction Hamiltonian Hrr as 

H S r = Y c i xx o- 

j 


Here, the reservoir is modeled as a bath of harmonic oscillators with Cj as the coupling constant. The dynamics of the system 
harmonic oscillator is obtained by tracing over the reservoir degrees of freedom. The optical tomogram from the Wigner function 
can be obtained using j l 1] 

u> (X, 9) = f W (X cos 9 — p sin 9, X sin 6 + p cos 9) dp , (43) 

where W ( x , y. t) is the Wigner function. Similarly, the corresponding Wigner function can also be reconstructed from the 
tomogram by inverse Radon transformation. The analytic expression of the tomogram for the system, initially in the coherent 
state |/3), is 


u(X,0,t) = 


y/(2NM+l)-{rMe~ 




(44) 


Here, k « is the Boltzmann constant and N = e xp(tm fc /fcBT)-i ' s t ^ lc average thermal photon number of the environment at 
temperature T. Also, r is the bath squeezing parameter, Re[u] denotes the real part of u and M = 1 — exp (—2 kt ), where k is 
the dissipation coefficient, analogous to the spontaneous emission term. 

In the corresponding figures for the tomogram with specific values of different parameters, we can observe the decay of the 
tomogram. Specifically, Fig. [9}a shows the tomogram of an initial coherent state, where we can see a beautiful valley like 
shape surrounded by a mountain. Interestingly, a similar tomogram has been observed for a binomial state of large dimension 
(cf. Fig. 2 in ll79ll ). However, in Fig. [9}b and c, we can see this sharp structure gradually fade away due to interaction with 
its surrounding. Thus, with increase in temperature, the effect of decoherence and dissipation, due to the ambient environment, 
deteriorates the obtained tomogram. Further, Fig. [TO] illustrates the effect of change of average thermal photon number and 
squeezing parameter, where in smooth (blue) and dot-dashed (cyan) lines we can observe the enhancement of decay. Similarly, 
dashed (red) and dotted (magenta) lines also show the effects due to changes in bath parameters for another set of parameters. 
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Figure 10: (Color online) The tomogram of a dissipative harmonic oscillator varying with time is shown for the initial coherent state parameter 
/3 = 2 and 9 = The smooth (blue) and dashed (red) lines correspond to the tomogram for N = 5, r = 1 for X = 1 and 2, respectively. 
Similarly, the dot-dashed (cyan) and dotted (magenta) lines correspond to the tomogram for N = 10, r = 4 for X = 1 and 2, respectively. 


VII. CONCLUSION 

Tomography is a powerful quantum state reconstruction tool. Its wide applicability in obtaining quasidistribution functions, 
quantum process tomography and density matrix reconstruction in quantum computation and communication is already estab¬ 
lished. However, these properties can get effected by the influence of the ambient environment. Here, an effort has been made 
to study the evolution of tomograms for different quantum systems, both finite and infinite dimensional, under general system- 
reservoir interactions, using the formalism of open quantum systems. The effect of the environment on the finite dimensional 
quantum systems, both spin and number-phase states, is to randomize the tomogram. For spin quantum states, single and two 
spin-i states are considered with open quantum system effects. For the number-phase states, a general expression is obtained 
and is illustrated through the example of a three level quantum (qutrit) system in a spontaneous emission channel. The increase 
in temperature tends to decohere the tomograms while squeezing is shown to be a useful quantum resource. Besides this, a 
tomogram for a spin-1 pure quantum state is also obtained. Further, the tomogram for an infinite dimensional system, the ubiq¬ 
uitous dissipative harmonic oscillator, is also studied. The results obtained here are expected to have an impact on issues related 
to quantum state reconstruction in quantum computation, communication and information processing. 
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Appendix 1 


The elements of the density matrix (l24l) are 


Pee (t) = e 2Tt Pee (0) ■ 

Pssit) = e -( r + r i2)tp ss (o) + (i _ e -( r - r i2)t) e _(r+ri2)t pee (0), 


Paa if 


= e-( r " 


Tl2)t paa (0) + (1 - e-( r + ri2 )‘) e-^-^Pee (0) , 


Pgg (t) = Pgg (°) + 0-~ e {r+r12) ‘) Pss (0) + (l - e 


-(r-ri 


+ 


(r+r 12 ; 

2r 


+ SSfi {(1 - «- <r - r “ ) *) - E# 2 (1 - =- 2r ')} 

Pes (t) = e -*(“0- ni 2)* e -3( 3r + ri2 ) t p es (0) , 

PeaW = e-^o+ ni2 ) t e-^ 3r - ri2 )Vea(0), 


(r-r 12 ) 


= -(r+ri 2 )t 


:) ‘) Paa (0) 


(1 - e -( r - r -) 4 5 6 7 ) 


(r-r 12 ) 

2 


Pee (0) ? 


if 

Psa f) — ^ 


_ g— 2iu)atg—Ft 
- p—2iQ,i2t /=) —Ft 


Peg (0) , 
e~ vt Psa ( 0 ) , 


p sg (t) = e _i (“ 0+ni2 )‘e _ 5( r + r i2) t 


Pag (0) + ({ 2f7l2e rt sin ( 2 °12 i ) 


+ r (l — e rt cos (2Sli 2 f))} + i { 2 fli 2 (l — e ri cos (2Qi 2 i)) — Te rt sin (2fli 2 i)} p es (0)] 
Pagit) = e-^o-n 12 )i e -i(r-r 12 )t p ag (0) - ( ff ({2f2 12 e~ r * sin (2 (l 12 t) 

+ T (l — e~ Tt cos (2f2i 2 f)) • — i {2S7i 2 (l — e~ Tt cos (2f2i 2 f)) — re _rt sin (2f2i 2 f)}) p ea (0) 
Here, all the matrix elements are written in the dressed state basis, which is connected with the bare state basis by 

Ip) = 1st) Is 2 >, 

l*> = ^(|ei)|P 2 ) + |pi)|e 2 )), 

|a) = ^(|e 1 )|fl 2 )-|ff 1 )|e 2 )), 

|e) = |ei)|e 2 >. 


Further, 




i - (p-fijY 


cos(k 0 rij) 

kgrij 


1 - 




cos(fc 0 rij) 
( k 0 r i:j ) 3 


(45) 


where jl = fi\ = /r 2 are the unit vectors along the atomic transition dipole moments, f,; 7 = f j — f ; , and fc 0 = ^ with 
wo = a)1 2 2 i the spontaneous emission rate is 

p _ U1 3 Ji 2 

1 i ~ 3t refic 3 ’ 

and the collective incoherent effect due to the dissipative multi-qubit interaction with the bath is 

Tij = Tji = ^/TiTjF (k 0 rij), 

for i ^ j with 


F{k 0 rij) = | l-(p,-fij) 


sin(konj) 

korij 


1-3 \ 


(korij) sin(korjj) \ 

(^orij ) 2 (korij ) 3 J 


Further, for the case of identical qubits, as considered here, fli 2 = fl 2 i, Ti 2 = r 2 i, and Ti = T 2 = T. 
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